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§ ■ Abstract 

, We describe a free-electron laser (FEL) in the Self Amplified Spontaneous Emission (SASE) 

regime quantizing the electron motion and taking into account propagation effects. We demonstrate 

' quantum purification of the SASE spectrum, i.e., in a properly defined quantum regime the spiking 

Q . behavior disappears and the SASE power spectrum becomes very narrow. 
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The Self Amplified Spontaneous emission (SASE) regime for a free-electron laser (FEL) is 
made up of three basic ingredients: high-gain, propagation or "slippage" effects and start-up 
from noise l|. The classical steady-state high-gain regime of FELs, with universal scaling 
and the introduction of the p-BPN parameter, was analysed in ref.|2], where the possibility 
of operating an FEL in the SASE regime was described. The first experimental observation 
of the high-gain regime, also starting from noise, was carried out in the microwave range 
using a waveguide in the Livermore experiment [3]- Presently, short wavelength FELs which 
amplify incoherent shot noise via SASE are of great interest worldwide as potential sources 
of ultra bright coherent X-ray radiation jj, 0] . 

Many theoretical studies of high-gain FELs do not take into account propagation 
effects and the initial noise is described by a small input signal or a small bunching. Other 
treatments assume that SASE is just steady-state amplification starting from noise, ignoring 
propagation effects 3, Ol • That approach does not give the correct temporal structure and 
spectrum of the SASE radiation as described in ref. 111. 

nn 

In ref. p, 110| it was shown that due to propagation there exists not only the steady-state 
instability of ref. Q], but also a superradiant instability, with peak intensity proportional 
to n^, where n is the electron density. This instability originates in the region near the 
rear edge of the electron bunch, producing a soliton-like pulse which grows and narrows 
as it slips over the electrons, preserving a hyperbolic secant shape with some ringing 
We stress that the mathematics and the physics of this superradiant instability, which is 
at the heart of SASE ^|, is completely different from the usual steady-state instability. A 
striking example of the difference between the steady-state and superradiant regimes is that 
superradiantly amplified spikes occur when the system is detuned from resonance, whereas 
the steady-state instability does not occur 0, Q]- In this case, the superradiant instability 
is actually stronger than for the case of exact resonance, since the radiation pulse propagates 
over fresh electrons, which have not been perturbed by the steady-state instability. 

As shown in ref. ij, a SASE FEL radiates a random series of spikes since, roughly 
speaking, the electron bunch contains many cooperation lengths Lc which radiate randomly 
and independently from one another. The final result is an almost chaotic temporal pulse 
structure with a broad spectral width. The number of spikes in the high-gain regime cor- 
responds approximately to the number of cooperation lengths in the electron bunch. Only 
when the length of the electron bunch is a few cooperation lengths is a partial "cleaning" 
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of the temporal profile of the radiation obtained, as described in ref. Q|. Furthermore, 
the total radiated energy does not behave as that predicted by steady-state theory. After 
reaching the steady-state value of saturation, the energy for SASE continues to increase due 
to the fact that the superradiant spikes do not saturate. However, SASE has one drawback 
with regard to its application as a useful source of short-wavelength coherent light : at 
short wavelengths many cooperation lengths Lc lie within the electron bunch. This implies 
a quasi-chaotic temporal structure of the radiation pulse and a consequent wide spectrum. 

Several schemes have been proposed in order to avoid the large spectral width associated 
with classical SASE. One of these involves a multiple wiggler scheme with a coherent seed 
laser 12, Q|. This method is based on the fact that the PEL interaction creates a strong 
bunching not only on the fundamental but also on the higher harmonics (even and odd), as 
shown analytically in 12, Hence, if the electron beam is injected in a second wiggler 
tuned to one of the higher harmonics, even or odd, the system starts radiating with intensity 
oc . This superradiant harmonic generation (SRHG) eventually evolves into the high-gain 
exponential regime. A difficulty in producing SRHG is that "it is necessary to make the 
first wiggler section long enough to bring the harmonic bunching well above noise, but not 
so long as to induce too much energy spread" |l^. This optimization of SRHG, termed 
high-gain harmonic generation (HGHG), has recently been performed by other authors jisl ]. 
However, this method presents difficulties when going to very short wavelengths, as pointed 
out in ref. due to the fact that each stage amplifies not only the coherent signal but 
also the noise. 

In this letter we propose a novel method for producing coherent short wavelength radia- 
tion with SASE. We introduce a quantum description of SASE which depends on a dimen- 
sionless quantum PEL parameter, p, which determines the number of photons per electron 
and/or the electron recoil in units of the photon momentum. We show that when p ^ 1 the 
SASE PEL behaves classically, i.e. in agreement with the SASE classical model. However, 
when p < 1 we obtain a quantum regime with features completely different from those of 
the classical regime, and which we shall refer as Quantum SASE. A surprising feature of this 
regime is the phenomenon of "quantum purification" , in which the random spiking behavior 
almost disappears and a strong narrowing of the spectrum occurs. Here we generalize a 
previous model 17] including propagation effects via a multiple scaling method used in the 
classical PEL theory of ref. [l^. This allows us to easily take into account the existence of 
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two different spatial scales in the FEL: the variation of the electron distribution on the scale 
of the radiation wavelength (describing the bunching) and the variation of the field envelope 
on the much longer scale of the cooperation length. 

The quantum FEL is described by the following equations for the dimensionless radiation 
amplitude A{z, zi) and the matter wave field \E'(^, z,Zi): 



^^^(^gl^ = -^^^^{e,z,z,)-zp[A{z,z^)e^'-c.c.]^{e,z,z,) (1) 

dA{z,z^) dA{z,zi) 1 ,,2 -ie n /n^ 

TT de\'^{e,z,zi)fe + t6A{z,zi). (2) 



dz dzi 27r Jq 

These equations were derived in detail in ref.|l7|, to describe Collective Recoil Lasing (CRL) 
in FELs and in atomic systems. In the CRL equations, the electrons are described by 
a Schrodinger equation for a matter-wave field \l/ in a self-consistent pendulum 

potential proportional to A, where | = |ap/ (Np), |ap is the average number of photons in 
the interaction volume V, and is the space-time dependent electron density, normalized 
to unity. In Egs. (ipi and Q we have adopted the universal scaling used in the classical 
FEL theory [ll. S lloj|. i.e. 6 = {k + kyj)z — ckt is the electron phase, where k^ = l-n jX^ and 
k = uj / c = 2tt/\ are the wiggler and radiation wavenumbers, z = z/Lg is the dimensionless 
wiggler length, Lg = Xw/Anp is the gain length, zi = {z — Vrt)/PrLc is the coordinate along 
the electron bunch moving at the resonant velocity Vr = c(3r = ck/{k + kw) and Lc = X/Airp is 
the cooperation length or coherence length, p = (a^ / Ack^Y^^ {e^n/ meoY^^ is the classical 
BPN g ^ V (A/2A.)(l + <) . the _ ene.gy .n ,„c= unit, a„. is the 

wiggler parameter and n is the electron density. Finally, p = (7—70) / P7o is the dimensionless 
electron momentum and 5 = (70 — 7r) / P7o is the detuning parameter, where 70 ~ 'jr is the 
initial electron energy in mc^ units. It can be easily shown that Eq. implies: 

^ j\e\^{e,zuz)f = 0. (3) 

Hence, the dimensionless density profile /o(^i) = Jq^ d6\'^\'^ is independent of z. This means 
that the spatial distribution of the particles does not change appreciably on the slow scale 
zi during the interaction with the radiation. 

Whereas the classical FEL equations in the above universal scaling do not contain any 
explicit parameter (see ref. the quantum FEL equations ((T)) and (j21) depend on the 



quantum FEL parameter 



From the definition of A, it follows that p\A\'^ = \a\'^/N is the average number of photons 
emitted per electron. Hence, since in the classical steady-state high-gain FEL A reaches a 
maximum value of the order of unity, p represents the maximum number of photons emitted 
per electron, and the classical regime occurs for p ^ 1. Note also that in Eq. (Q) p appears 
as a "mass" term, so one expects a classical limit when the mass is large. As we shall see, 
when p < 1 the dynamical behavior of the system changes substantially from a classical to 
a quantum regime. 

A further inspection of Eq. shows that \1/ depends explicitly on 6, which describes 
the variation (bunching) on the scale of the radiation wavelength A. The parametrical 
dependence on zi is induced by the slow spatio-temporal evolution of the field amplitude 
A. This evolution is described by Eq. (j2)), in which the bunching factor is the ensemble 
average of e~*^. In other words, the A^-particle system is described as a quantum ensemble 
represented by the matter- wave field \E'. This model has been used previously to describe the 
quantum regime of an FEL and of the Collective Atomic Recoil Laser (CARL) 
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2J,|2^|2^, neglecting the dependence on zi, i.e., propagation. This is appropriate for the 
FEL when slippage due to the difference between the light and electron velocities is neglected, 
which is never true in the SASE regime where the propagation from one cooperation length 
to another is fundamental. 

The classical limit of the FEL can be explicitly shown as follows. Eq.(P]) can be trans- 
formed into an equation for the Wigner function, as shown in ref. 

dW{e p,z,,z) dWie,p z,,z) ^^^.e ^ ^.^-.e. 
oz oO ^ ' 



^ (5) 



whereas Eq.(j2} becomes 

dA dA 1 '•+°° '•2- 

dz 



dA 1 /""^°° f'^'^ 
+ 7^ = —/ dp deW{e,p,zi,z)e-'' + t6A. (6) 
dzi 2n Jq 

In the right hand side of Eq. (jSJ, the incremental ratio p[W{9,p + l/2p, zi, z) — W{9,p — 
l/2p,zi,z)] — > dW{9,p, zi, z)/dp when p — > oo. Hence, for large values of p, Eq. (0) 
becomes the Vlasov equation: 

dW{9,p, zi,z) ^ _ dW{9,p, zi,z) _ ^ A*e-'^) ^^(^'^' ^i'^) = q (7) 
dz d9 ^ ^ dp ' 

Eqs. (0) and © provide the description of the CRL model in terms of the Wigner function 

for the electrons, whereas Eqs. (jH)) and ((Tj) are equivalent to the classical FEL model of 
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ref. Note that Eqs. (jH)) and ((Tj) do not depend explicitly on p, as must be the case in 



the classical model with universal scaling 




10|. We briefly mention that Eq. (0) for the 



Wigner function has a broader validity than the Schrodinger equation (^Q), since it can also 
describe a statistical mixture of states which cannot be represented by a wave function but 
rather by a density operator. 

Eqs. and (j21) are conveniently solved in the momentum representation. Assuming that 
"^{6, zi, z) is a periodic function of 6', it can be expanded in a Fourier series: 

oo 

m{e,z,,z)= J2 Cn(^i,^)e*"' (8) 

n=— oo 

SO inserting Eq. (jH)) into Eqs. (H)) and we obtain 

= -^-^Cn- p{ACn-l- A*Cn+l) (9) 

OA OA ^ ^ - ^ 

n=— oo 

Eqs. and (jlUj) are our working equations and their numerical analysis will be discussed in 
the following. Note that from Eq. (jS)) it follows that |c„p is the probability that an electron 

has a dimensionless momentum p = n/p (i.e. mc(7 — 70) = {hk)n), 

00 

b = (11) 

n=— 00 

is the quantum expression of the bunching parameter and 

^ 00 

{P) = -_Y. (12) 
P ^ 

n=— 00 

is the average dimensionless momentum. Note that the quantum bunching (jlip requires a 
coherent superposition of different momentum states. The stability analysis of Eqs. © and 



(I10|) has been carried out in ref. 2J, |25|, |26[ for the case with no propagation. We assume 
that the system is in an equilibrium state with no fleld, A = 0, and all the electrons are in 
the state n = 0, with cq = 1 and Cm = for all m 7^ 0. This equilibrium state is unstable 
when the dispersion relation 

(A-5)(a^--^)+1 = 0, (13) 

has complex roots. Notice that this dispersion relation coincides with that of a classical 
FEL with an initial energy spread with a square distribution and width 1/p [l^, i.e., this 



extra term represents the intrinsic quantum momentum spread which, in dimensional units, 
becomes hk. In fig^we plot the imaginary part of A as a function of 6 for different values of 
p. The classical limit is obtained for p ^ 1 (see figCJi). In this case, the system is unstable 
for 5 < 2, with maximum instability rate ImA = at 5 = 0. When p is smaller than 

unity (figHh-f), the instability rate decreases and the peak of Im(A) moves around 6 = l/2p 
(i.e. mc^'jQ — 7r) = hk/2) with peak value ImA ==_v/p and full width on S equal to 4p^/^ 
(i.e. mcA'j = 4{hk)p^^'^). As discussed in ref. j^Sl' ^'^^ p ^ 1 the electrons have almost 
the same probability of transition from the momentum state n = Oton = lorri = — 1 
(i.e. |cip ~ |c_ip), absorbing or emitting a photon. On the contrary, in the case p < 1, 
|cip <^ |c_ip, i.e. the particles can only emit a photon without absorption. 

As shown by the linear stability analysis discussed above, the quantum regime occurs for 
small value of p, when an electron emits only a single photon. In this limit, the dynamics 
of the interaction is that of a system with only two momentum states, i.e. the initially 
occupied state with n = and the recoiling state with n = —1. In this limit, Eq. (jHl) and 
f|Tn|l . after defining the "polarization" S = coc*_iexp[i6z] and the "population difference" 
D = |cop — |c_ip, reduce to the so-called "Maxwell-Bloch equations" for a two-state system 

Eh 

d 

-g^S{zi,z) = -iAS{zi,z) + pA{zi,z)D{zi,z) (14) 

^D{zu z) = -2p [A(zi, zYS{z,, z) + c.c] (15) 

^ + ^ = ^(-1,4 (16) 

where A = 5 — l/2p and A = Aexp {—i6z). Notice that the parameter p in Eqs. ()14|) - 
(Uni) may be eliminated by redefining the variables as A' = y/pA, z' = ^z., z[ = y/pzi 
and A' = A/^/p. With this quantum universal scaling the cooperation length becomes 
L'^ = Lc/ y/p oc 1/ ^/n. An interesting result of this scaling is that the CRL model of Eqs. 

and can now be interpreted as a Schrodinger equation for a single particle with a 
"mass" p^/^ in a self-consistent pendulum potential. This provides an intuitive interpretation 
of the classical limit which holds when the particle's "mass" is large. 

We now discuss some analytical results arising from the quantum SASE model. It has 
been shown in a previous work that the classical model of Eqs. and (0) admits a 
soliton-like self-similar solution A{z,Zi) = ZiAi{y) when 6 = 0, where y = yfz{{z — Z]) ^ 
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{z — VrtY^'^{ct — z)/L^J'^ (where we have assumed /5r ~ 1) This means that the radiation 
pulse propagates over the electron bunch (i.e. at different values of Zi), preserving its shape 
but increasing with amplitude oc zi and narrowing with width oc l/^/zi. It is possible to 
show that in the linear regime the radiation pulse has a maximum at zi = z/3 and 
propagates at a constant velocity Vg = 3vr/{2 + [3r). 

An analogous self-similar solution also exists for the quantum equations (jl4|) - (jl6|) when 
A = [281, i.e. A{z,zi) = ziA2{x), where x = pzi{z — zi) = {z — Vrt){ct — z)/L'^ (with 
L'^ = Lcl\U)). The shape of the radiation pulse is similar in the classical and quantum 
cases, but in the quantum case its width decreases as 1/pzi. In this case, the radiation 
pulse in the linear regime has a maximum ai zi = z/2 and moves at the constant velocity 

From the features of the self-similar solutions we can deduce some important aspects of 
the nature of the random spikes emitted in the SASE regime. First, both in the classical 
and in the quantum regime the dimensionless amplitude A of the field is proportional to 
Zi (X p (X A^^/^, so that the number of emitted photons |ap = A^p|y4p oc A^^, ie. is 
superradiant. Whereas the characteristic spatial length in the classical regime is Lc, in the 
quantum regime the characteristic length is L'^ = Lc/ y/p ^ Lc for p ^ 1. 

Let us now briefly restate the reasons for the classical random spiking behavior. If the 
bunch length contains many cooperation lengths L^, each of them radiates a superradiant 
spike independently as the electron bunch propagates into the wiggler. Each spike is only 
roughly represented by a self-similar solution because the radiation from one portion does 
not find fresh, electrons but electrons which have already interacted with radiation. As 
discussed in |1| the number of spikes increases with the number of cooperation length in the 
bunch, i.e. Lh/2nLc . 

We now describe the numerical solution of Eqs. and (j21) which demonstrates the differ- 
ent dynamical behavior in the classical and quantum SASE regimes. The SASE simulation 
has been performed assuming all the electrons are initially in the n = momentum state 
and that there is a weak randomly fluctuating modulation in the electron density along the 
electron pulse, as is appropriate to model random electron shot noise. The initial condi- 
tions for all the simulations are therefore A{zi,z = 0) = 0, c^i(zi,z = 0) = b^e^'^^^'^^ and 



Co{zi,z = 0) = a/1 — &0) where 60 = 0.01 and (f){zi) is a randomly fluctuating phase with 
values in the range [0, 27r). Since there is not a radiation seed, we assume 6 = 0. 
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Fig. 2 shows the numerical solution for Lh = 30Lc and ^ = 50 for the classical regime 
(p = 5, left column) and the quantum regime (p = 0.05, right column). Figures |2K and|2b 
show the intensity |Ap as a function of the dimensionless variable Zi (i.e. the coordinate 
along the electron bunch in units of Lc in the electron rest frame), whereas figures |2t and 
|21i show the power spectrum P{uj) as a function of = {Lc/c){uj' — u), where u' is the field 
frequency and is the carrier resonant frequency. In fig. andl^D , = is the trailing 
edge and zi = 30 is the leading edge of the electron bunch. Therefore, the region on the left 
of the dotted line is the radiation on the electron beam, on the right is free propagation in 
vacuum. 

Alternatively one can interpret fig. and|2b as the temporal behavior of the intensity 
for an observer which is at a given position in the wiggler and which will see the pulse as it 
appears from right to left in the figure. We stress that in the classical theory with universal 
scaling the two cases should be identical. On the contrary, the dramatic difference is evident 
from fig. 121 the temporal structure in the classical limit (fig. 12^) is almost chaotic and the 
width of the spectrum is large (fig. |2t)- Conversely, the temporal behavior in the quantum 
limit (fig. ISb) shows a purification of the initially noisy evolution so that the temporal 
structure looks similar to the self similar solution one would obtain with a coherent seed 
signal. The corresponding spectrum becomes extremely narrow (fig. |21i), much sharper than 
the classical one. It can be seen in fig. |211 that the frequency is shifted by l/2p = 10, 
in agreement with the predictions of the linear analysis described earlier (see fig. The 
behaviour of the system is similar to what would be expected if the quantum cooperation 
length is much larger than the classical cooperation length, so that all electrons radiate 
coherently in the quantum regime. The difference between classical and quantum SASE 
behavior is confirmed by a comparison of fig. and Eb, which show the total radiated 
energy, E (in units of |Ap), in the classical and quantum regimes respectively, where 

rz+L 

E= \A\'^dzi. 
Jo 

In the quantum regime we observe a behavior similar to that which one would obtain with 
a coherent seed signal in the long bunch superradiant case 0, Q . 

Finally, we discuss the reason for quantum purification of the spectrum. As remarked 
earlier, in fig ^ the gain bandwidth (which is the reciprocal of the "real" coherence length) 
decreases as p^^^ in the quantum regime. Note that the cooperation length in the quantum 
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regime, L'^ — Ld^i has the same dependence on p. Hence, one can understand that in 
quantum SASE, i.e., for small value of p, the system behaves as if the startup of the PEL 
interaction is initiated by a coherent bunching or coherent seed. 

In conclusion, in this letter we have given a proof of principle of the novel regime of 
Quantum SASE, with dynamical properties very different from "normal" classical SASE. 
In particular, quantum SASE predicts quantum purification of the temporal structure and 
spectrum. The possibility of experimental observation of this quantum regime is under 
investigation and will be discussed elsewhere. 

We acknowledge useful discussions with S. Bertolucci, L. Serafini, M. Ferrario and L. 
Palumbo. This work is supported by INFN. 
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FIG. 1: Imaginary part of the unstable root of the cubic equation (fT^ vs. 6, for 1/2/3 = 0, (a), 
0.5, (b), 3, (c), 5, (d), 7, (e) and 10, (f). 
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FIG. 2: Numerical solutions of eq. Q and in the classical and quantum regimes : Graphs (a) 
and (b) show the scaled intensity |^p(zi) at z = 50 when the system evolves classically (p = 5) 
and quantum mechanically (p = 0.05) respectively. Graphs (c) and (d) show the corresponding 
scaled power spectra P{0) at z = 50 when p = 5 and p = 0.05 respectively. The frequency shift in 
(d) is l/2p = 10 in agreement with that predicted from fig. ^ In all cases, Z = 30 and 6 = has 
been used. 
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FIG. 3: Scaled energy, as a function of z when (a) p = 5 and (b) p = 0.05. In both cases 6 = 
and L = 30. 
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